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INTRODUCTION 

The production of heavy clusters in nuclear reactions is 
important in a wide variety of applications: radiation shielding, 
space engineering and design, medical physics, accelerator 
design, and more. According to the Cascade Exciton Model 
(cem) mm , there are three ways high-energy heavy clus¬ 
ters can be produced. The first way is via coalescence of 
nucleons produced in the IntraNuclear Cascade (INC). The 
second way is via the preequilibrium model. The last way 
is via Fermi breakup. Previous work in CEM examines the 
impact of expansions of the preequilibrium model and Fermi 
breakup model on heavy cluster production [3-7]. The present 
work studies the impact of expanding the coalescence model 
on heavy cluster spectra. CEM03.03 Q , the default event gen¬ 
erator in the Monte Carlo N-Particle transport code version 6 
(MCNP6) OD for intermediate-energy nuclear reactions, is ca¬ 
pable of producing light fragments up to 4 He in its coalescence 
model. In the present study, we have expanded the coalescence 
model to be able to produce up to 7 Be. Preliminary results are 
promising. 

BACKGROUND 

When the cascade stage of a reaction is completed, CEM 
uses the coalescence model described in Ref. ora to “cre¬ 
ate” high-energy d, t, 3 He, and 4 He by final-state interactions 
among emitted cascade nucleons outside of the target nucleus. 
The coalescence model used in CEM is similar to other coa¬ 
lescence models for heavy-ion-induced reactions. The main 
difference is that instead of complex-particle spectra being esti¬ 
mated simply by convolving the measured or calculated inclu¬ 
sive spectra of nucleons with corresponding fitted coefficients, 
CEM03.03 uses in its simulations of particle coalescence real 
information about all emitted cascade nucleons and does not 
use integrated spectra. (Note that the coalescence introduced 
recently in the Liege intranuclear cascade (INCL) [11-14], is 
in a way similar to the coalescence considered by CEM as 
proposed in Ref. mm, with the main contrast being that 
INCL considers coalescence of INC nucleons on the border of 
a nucleus, inside the target-nucleus, while CEM coalesces INC 
nucleons outside the nucleus.) We assume that all the cascade 
nucleons having differences in their momenta smaller than p c 
and the correct isotopic content form an appropriate composite 
particle. The coalescence radii p c as used in CEM03.03 are: 
Incident energy, T, < 300 MeV or > 1000 MeV 

p c (d) = 90 MeV/c ; 

Pc(t) = Pc( 3 He) = 108 MeV/c ; (1) 

p c { A He) = 115 MeV/c. 


300 MeV < T < 1000 MeV 

p c (d) = 150 MeV/c; 

p c (t) = p c ( 3 He) = 175 MeV/c ; (2) 

p c ( 4 He) = 175 MeV/c . 

If several cascade nucleons are chosen to coalesce into 
composite particles, they are removed from the distributions 
of nucleons and do not contribute further to such nucleon 
characteristics as spectra, multiplicities, etc. 

COALESCENCE EXPANSION 

The Coalescence Model in CEM03.03 allows for coales¬ 
cence up to 4 He. We have expanded this to additionally allow 
for the coalescence of 6 He, 6 Li, 7 Li, and 7 Be. 

CEM03.03 uses the simplest version of the coalescence 
model ora and checks only the momenta of nucleons emit¬ 
ted during the cascade stage of reactions, without checking 
their coordinates. 

The momentum, p , of each nucleon is calculated relativis¬ 
tic ally from its kinetic energy, T (CEM03.03 provides in its 
output files the energy of particles, but not their momenta). 
Coalescence occurs if each nucleon in the group has |Ap| < p c , 
where Ap is defined as the difference between the nucleon 
momentum and the average momentum of all nucleons in the 
group. 

The coalescence model of CEM03.03 first checks all nu¬ 
cleons to form 2-nucleon pairs, their momenta permitting. It 
then checks if an alpha particle can be formed from two 2- 
nucleon pairs (either from 2 n-p pairs or from a n-n and p-p 
pair). After this it checks to see if any of the 2-nucleon pairs 
left can combine with another nucleon to form either tritium 
or 3 He. And lastly, it checks to see if any of these 3-nucleon 
groups (tritium or 3 He) can coalesce with another nucleon to 
form 4 He. 

The expanded coalescence model then takes these 2- 
nucleon pairs, 3-nucleon (tritium or 3 He only) groups, and 
4 He to see if they can coalesce to form heavy clusters. 4 He 
can coalesce with a 3-nucleon group to form either 7 Be or 
7 Li. Two 3-nucleon groups can coalesce to form either 6 Li 
or 6 He. And 4 He can coalesce with a 2-nucleon pair to form 
either 6 Li or 6 He. All coalesced nucleons are removed from 
the distributions of nucleons so that our coalescence model 
conserves both atomic- and mass-numbers. In addition, this ex¬ 
pansion requires an insignificant amount (2-3 %) of increased 
computation time. 

For additional details of the Coalescence Model expan¬ 
sion, see Ref. E3. 


Coalescence Parameter p c 

As mentioned above, p c determines how dissimilar the 
momenta of nucleons can be and still coalesce. p c was ex¬ 
panded to also include a value for heavy clusters, or light 
fragments (LF): p c (LF). Our new p c ’ s for incident energies, 
T, less than 300 MeV or greater than 1000 MeV are: 

p c (d) = 90 MeV/c ; 

p c {f) = pcfHe) = 108 MeV/c ; (3) 

PcC He) = 130 MeV/c. 

p c (LF ) = 175 MeV/c . 

And for 300 MeV < T < 1000 MeV 

p c (d) = 150 MeV/c; 

p c (t) = p c CHe) = 175 MeV/c ; (4) 

Pci 4 He) = 205 MeV/c . 

Pc(LF) = 250 MeV/c . 

Note that the pd 4 He) was also increased compared to 
the old p c values. Too many alpha particles were lost (coa¬ 
lesced into heavy clusters), and therefore we compensated by 
increasing the coalescence of 4 He. 

RESULTS AND ANALYSIS 

Examples of results of our coalescence expansion are 
displayed in Figs. 1-3. The upgraded CEM03.03F without 
coalescence expansion (blue solid lines) and the upgraded 
CEM03.03F with coalescence expansion (red dashed lines) 
are compared with experimental data (green circles). This 
coalescence expansion analysis is part of a larger project [3-7] 
aimed at producing high-energy light fragments in spallation 
reactions. CEM03.03F refers to the upgraded CEM03.03 code, 
which has been upgraded with a Modihed-Exciton-Model ex¬ 
pansion and a total reaction cross section model improvement. 
The blue solid lines contain both of these improvements over 
the original CEM03.03, and the red dashed lines contain both 
of these improvements plus the coalescence expansion. 

Fig. [I] displays fragment production spectra of 6 He and 
6 Li for the reaction 1200 MeV p + 197 Au. Experimental data 
by Budzanowski et al. m (green circles) are compared with 
results from CEM03.03F without coalescence expansion (blue 
solid lines) and CEM03.03F with coalescence expansion (red 
dashed lines). The coalescence expansion increases the pro¬ 
duction of high-energy 6 He and 6 Li, and improves agreement 
with experimental data. 

Fig. [5] displays the fragment production spectra of 7 Li 
and 7 Be for the reaction 480 MeV p + nat A g. Experimental 
data by Green et al. ini (green circles) are compared with 
results from CEM03.03F without coalescence expansion (blue 
solid lines) and CEM03.03F with coalescence expansion (red 
dashed lines). Again, the coalescence expansion increases the 
production of heavy clusters, and improves agreement with 
experimental data. 

These reactions also highlight how the coalescence can 
produce heavy clusters not just of high-energy, but also of 


low- and moderate-energy, thus improving agreement with 
experimental data in these energy regions as well. 

Fig. [3] displays experimental results of the reaction 480 
MeV p + naT Ag —> 6 Li by Green et al. fTTI (green circles), 
compared with simulations from results from CEM03.03F 
without coalescence expansion (blue solid lines), CEM03.03F 
with coalescence expansion (red dashed lines), and the orig¬ 
inal CEM03.03 (brown dashed-dotted lines). Even without 
the coalescence expansion, CEM03.03F (which contains a 
preequilibrium expansion and a total reaction cross section im¬ 
provement) yields much better results than CEM03.03 without 
any of these improvements. Adding the coalescence expansion 
produces even better results. 

Similar results for many other reactions induced by pro¬ 
tons, neutrons, and heavy ions (the last are simulated with 
LAQGSM03.03 0, but with an extended coalescence model 
as described in Ref. lfl8l ) and further discussions can be find 
in Refs. lfl5lH8l . 




Fig. 1. Comparison of experimental data by Budzanowski, et 
al. El (green circles) for the production of 6 He and 6 Li at 
an angle of 20° from the reaction 1200 MeV p + 197 Au, with 
results from CEM03.03F without coalescence expansion (blue 
solid lines) and CEM03.03F with coalescence expansion (red 
dashed lines). 






















Fig. 2. Comparison of experimental data by Green et al. fT71 
(green circles) for the production of 7 Li and 7 Be at an an¬ 
gle of 40° from the reaction 480 MeV p + nat Ag , with re¬ 
sults from CEM03.03F without coalescence expansion (blue 
solid lines) and CEM03.03F with coalescence expansion (red 
dashed lines). 

CONCLUSIONS 

Expanding the coalescence model within CEM yields 
increased production of heavy clusters in nuclear spallation 
reactions, particularly in the high-energy region, but also in 
the low- and moderate-energy regions. Preliminary results in¬ 
dicate this coalescence expansion yields improved agreement 
with experimental data. We plan to implement these upgrades 
into MCNP6. 
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